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The entropy production occurring in relativistic hydrodynamical systems 
such as the quark-gluon plasma (QGP) formed in high-energy nuclear colli- 
sions is explored. We study mechanisms which change the composition of the 
fluid, i.e. particle production and/or chemical reactions, along with chemo- 
and thermo-diffusion. These effects complement the conventional dissipative 
effects of shear viscosity, bulk viscosity, and heat conductivity. 



Introduction. There is fundamental interest in the study of high energy density matter, 
which promises to illucidate further the properties of the strongly interacting vacuum and 
the high temperature phases of strongly interacting matter |jl], |^]. The method of choice 
for the description of the rapidly evolving and possibly inhomogeneous quark-gluon plasma 
(QGP) has been often the hydro dynamical approach |3|, |^. It is generally believed that 
QGP existed shortly after the Big Bang initial phase of the universe, and that it may be 
re-created in relativistic nuclear collisions ||5[ |6|. 

In order to gain a better understanding of this new form of matter, we explore here 
the entropy production due to numerous dissipative effects arising in the evolution of QGP 
towards the point of hadronization. The entropy content of the QGP phase at hadronization 
determines the observable particle multiplicity 0]. Therefore, in order to establish via 
particle multiplicity the initial conditions reached in heavy-ion collisions the effect of 
entropy production in QGP evolution has to be understood. 

To begin with, we present here the main result of this paper, before embarking on its 
derivation. The entropy production is given by: 

i id 

-^Q.A^^'^Qu + |>Vf A^W^ A? + Cid.u^f , (1) 

see Eq. (|2^). We will show that the four-divergence of the entropy current, d^s^, is always 
nonnegative and vanishes only, when all space-time gradients of temperature T, of chemical 
potentials /Xj, and of four- velocity vanish, and particle densities of species "i" assume 
their chemical equilibrium values. - The projector A is: A^^, = g^^ — Ufj_u,^, with g^^ = 
diag(l, —1, —1, —1) and u'^ = 1; we use units such that fi = c = ks = 1- 

The first two terms on the right-hand side of Eq. (|^) are new. They are due to the sources 
i7j of particle production and chemo- and thermo-diffusion contributions involving gradients 
of chemical potentials and temperature , respectively. These terms arise in the continuity 
Eq. (^) below; a denotes the matrix of mutual diffusion coefficients between the chemical 
QGP components. Furthermore, the constants k, t], and C denote the transport coefficients 
of heat conductivity, shear viscosity, and bulk viscosity, respectively. The corresponding 
terms in Eq. (|l|) involve the heat-flow four- vector Q^, Eq. 
Eq. (|l8|) , and have been obtained earlier ^ . 

Generally, dissipative processes are important, when space-time gradients of hydro- or 
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17), and the shear tensor W, 



thermodynamical quantities in the system become large relative to its relaxation scales. 
The latter and especially the transport coefficients and source terms introduced in Eq. (|^) 
have to be fit to experiment or calculated from an underlying microscopic transport theory 
p. For a recent extensive review of such state-of-the-art microscopic calculations and 
new gauge theory results we refer to Ref. [^] and its numerous references. 

In the QGP case, particle production converting kinetic energy of a nuclear collision 
into high multiplicity of new particles is a characteristic feature, together with annihilation 
processes, changing the composition of the plasma |11, 13|. Motivated by this fact, our 



aim presently has been to incorporate these dissipative and entropy producing mechanisms 
consistently into the relativistic hydro dynamical framework. 

Equations of Motion. Presently we consider the relativistic energy-momentum tensor 
for an imperfect fluid: 

T'^^ = (e + P)u^'u'' - Pg^''' + dT^'^ = T^^^ , (2) 

where e and P denote its local energy density and pressure, respectively. The additional 
term 5T^^^ incorporates the dissipative effects within the fluid @, ^. It will be obtained in 
detail in the following. We point out that the corresponding transport coefficients have to 
be determined consistently, i.e., including also all contributions from particle production 
or other composition changing processes, if present in the system, besides the scattering 



contributions P, IC]. 



The evolution of the hydrodynamical system then is governed by the equation of motion, 

d.T'^" = r , (3) 

where describes the density of external forces, gravity for example, which act on the fluid 
locally. A closed system has = 0, which results in total four-momentum conservation. 
In order to obtain a complete description of the fluid, the equation of motion (2) has to be 
supplemented by an equation of state relating the energy density to the pressure. In the 
simplest case of a noninteracting ultrarelativistic system, we have P = e/3, while a realistic 
QGP equation of state is given in Ref. [0]. 

Often, particularly if there exist conserved charges or particle numbers, a simple func- 
tional relationship between pressure and energy density does not exist. Then, assuming local 
thermal and possibly chemical equilibrium, both can still be represented parametrically in 



2 



terms of a temperature T = T{x) and possibly several chemical potentials /Uj = ^i{x). In 
this case, one has to complement Eqs. (|2|)-(^) by continuity equations, 

d^p>i = d^p^u^' = Q , (4) 

where pi = Pi{x) denotes the particle number density for species "i" in the local rest frame of 
a fluid cell. We will not distinguish the cases of charge or particle number conservation and 
count each (anti-)particle species separately. - The local rest frame is always determined 
by = u^Q^ = (1,0,0,0). Thus, by definition, we have p^^^^^ = pi. 

Similarly, the total energy density is defined to be equal to e. Therefore, we must require 
that dissipative terms do not contribute to T^^ in the local rest frame, 

Tjjj = e 6Tl>^^ = Uf.uJT'^'' = , (5) 

cf. Eq. (^). This constraint plays an important role in determining the general structure of 
ST^"^ consistently with the Second Law of Thermodynamics. 

In order to incorporate diffusion and particle production and/or chemical reactions 
among the constituents, which change the composition of the fluid, we now generalize 
Eq. (I) to read: 

= + A^'^a.7^i) = , (6) 

with the new four-current pi having the property u^p^ = pi, since u^A^'^ = 0. Furthermore, 
J'i denotes the local source for particle species "i", which is a functional of the densities, 
J'i{x) = J'i[pj]{x). It can be calculated in terms of microscopic reaction cross sections. 
We do not consider here the possibility of external sources. - In a homogeneous system 
such continuity equations become familiar chemical rate equations, which have previously 
been applied to strongly interacting relativistic systems on the hadronic and partonic level 
B 0,1161. 



The explicit form of the diffusion term oc dTZi in Eq. (|6|) will be determined consistently 
with the Second Law of Thermodynamics. - One might try to follow the nonrelativistic 
ansatz which recovers the usual diffusive contribution to the particle current in Eq. (^) in 
the local rest frame: 

T^i = Yl ^ijPj ' 3i ^ -VTZi = - ^ ^ij'^Pj , (7) 

j j 



3 



where dij represents the mutual diffusion coefficient between species "z" and "j". How- 
ever, this ansatz turns out to be incompatible with the requirement that the entropy never 
decreases and will be suitably generalized below. 

Entropy Production in an Imperfect Fluid. Here we formulate a consistent theoretical 
framework to study the entropy production due to dissipative effects within a relativistic 
fluid. We assume local thermal equilibrium and characterize the particle densities by ap- 
propriate chemical potentials. 

We remark that introducing chemical potentials does not necessarily imply chemical 



equilibrium in the sense of saturating the available phase space |17]. - For example, con- 
sider a thermalized electron-positron plasma with one chemical potential regulating the 
charge density, as usual, and with a number of e'''e~-pairs which does not correspond to the 
given temperature. Such a situation can be described by introducing a phase space occu- 
pancy factor in addition to the charge chemical potential or, equivalently, by two chemical 
potentials regulating the number of electrons and positrons separately. - We presently ad- 
mit such cases, recall Eq. (P), which typically involve charge conservation side by side with 
dynamically changing particle pair populations in the fluid. 

In order to obtain a relationship between the thermodynamical quantities and variables 
discussed so far and the entropy, in particular, we employ the (local) equilibrium relations, 

-PV = Q{T,V,^ii) = U-TS-J2l^^Ni , (8) 

i 

where U, S, and Ni denote the internal energy, entropy, and particle numbers, respectively. 

Two comments are in order here. Firstly, having nuclear collisions in mind, we employ 
the grandcanonical description despite the fact that the fireball created in a single event is 
thermodynamically closed. There is no heat or particle bath. However, instead of generating 
the grandcanonical ensemble by our system being connected to such a bath, we consider a 
large number of events to form the ensemble. The events have to be identical with respect to 
the conserved total energy-momentum, angular momentum, and baryon number, to name 
the most relevant "macroscopic" observables. Nevertheless, microscopically, the initial state 
may differ from event to event. Then, under the assumption that the entropy is maximized. 



the grandcanonical ensemble is obtained |17, 18|. In the hydrodynamical approach we apply 
this reasoning locally to each fluid cell. 

Secondly, consider finite size effects present in small volume systems. In this case, the 
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usual (y oo) statistical description breaks down, if the unphysical fluctuations of exactly 
constrained quantities are not eliminated. [| This issue needs to be addressed in our context 
in the future. 

We also need the Duhem-Gibbs relation for the densities corresponding to Eq. (^), 

Ts = e + P-Y,f^tPi ■ (9) 

i 

Then, using Eqs. and (^), we obtain by contraction with Ui, of the equation of motion 
(I): 

i 

i i 

= ^|,[T~S^ + Y,P^P'^]-U^^^P + Uu^^ST^^' , (10) 

i 

where we introduced an entropy four-current = su'' and the currents pf = piU^, as 
before. For the last equality, we made use of = 1. 

For conservative external forces which can be derived from a potential (density), = 
d^4), using Eq. (P) and n^9^ = d/dr, we can further simplify Eq. (|l^) [A = dA/dr): 

= T^^~s^' + Y,^l^{Ji-^^^^''^,'R{) + sf + Y.PiP^-P-'t> + Uu^^,5T^''' 

i i 

= Td^sf' + J2M-d,A^"dM+u,d^6T^''' , (11) 

i 

where the terms involving d/dr cancel identically. In order to demonstrate this cancellation, 
we multiply them by ^ • and employ Eq. (|8|), 

SdT + J2 Nidpi - VdP - Vdcp 

i 

= d{TS) + d{J2 PiNi) - d{PV) - TdS - ^ PidNi + PdV - Vd<j) 

i i 

= dU - TdS - E PidNi + PdV - Vdcj) = . (12) 

i 

The last equality expresses the First Law of Thermodynamics, which is incorporated here. 

^Such constraints arise due to exact conservation laws, unbroken internal symmetries, and quantization of 
the motion of the constituents (see, e.g., Refs. |]l9| , |2^, §2[). For example, the thermodynamical properties 
of small quark-gluon plasma droplets start to deviate strongly from the asymptotic {V oo) behavior, if 
the dimensionless parameter TV^^^ becomes of order one or less. 



5 



Finally, we introduce the proper entropy four-current defined by: 

= s^' +T-^u^ST^"' + Y,i^^^^"'dyni . (13) 

i 

This form is found by 'trial and error' and, by the constraint (^) and Uf^A^'^ = 0, implies 
■^("o) = •s in the local rest frame. 

We remark that the difficulty to identify a viable form of the entropy four-current, i.e. 
the admissible structure of 5T^^ , Ci, and TZi, resides in the requirements that the entropy 
density be equal to s in the local rest frame and that it never decreases anywhere in the 
system. In the following we will determine the function Ci consistently together with TZi 
and 6T'^'^ , in order to obtain the new contributions due to diffusion and, particularly, due 
to particle production, which have not been considered before. 

We combine Eqs. (0) and (|l3|), which yields: 

Td^s'' = -Y,f,iJi-T-\d^T)u,5T>''' + {d^u,)6T^'' 

i 

+ J2{f^i + TCi)d^A^'''d,n., + Tj2{d^Q)^^''d,'JZ^ . (14) 

i i 

This equation describes the entropy production in an imperfect fluid evolving under the in- 
fluence of conservative external forces, which do not generate entropy (Liouville's theorem). 
In order to agree with the Second Law of Thermodynamics for the closed system (/'^ = 0), 
the right-hand side of Eq. ( |l^ has to be nonnegative for all fluid configurations, such that 
d^s^ > always. This requirement is implemented in the following steps. 
First, we put a constraint on the sources introduced in Eq. (^), 

^^iJi < . (15) 

i 

This constraint and its implications will be further discussed below, beginning with Eq. (|2^). 
Under fairly general assumptions about the particle producing processes or chemical reac- 
tions the inequality ( p^ ) guarantees that the system is driven towards chemical equilibrium. 

Second, the structure of 5T^^ including terms of first order space-time gradients, which 
is subject to the constraint (^) and compatible with 5^s^ > 0, is well known [^: 

5T^''' = k{A''^u'' + A'^'^u'') + 7?A^^ A^^W^^ + (^'""d^u^ , (16) 

where k, ry, C > denote the coefficients of heat conductivity, shear viscosity, and bulk 
viscosity, respectively, and with the heat-fiow four-vector and the shear tensor W^u 
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defined by: 



= d^T-Tu'^d.u^ , (17) 
2 

y^i^u = d^Uy + dyu^ - -g^iud^u' . (18) 



Third, as we shall see in Eq. (23), the contribution to the right-hand side of Eq. (14) of 
the diffusion terms involving TZ becomes a nonnegative quadratic form if we identify: 

3 

where the symmetric matrix a consists of the mutual diffusion constants and supposedly 
has only nonnegative eigenvalues. We remark that IZi as defined in (|l9|) yields a diffusive 
particle current in the local rest frame, 

I ^ -V7^. = - ^ ^(V^, - ^VT) . (20) 

i 

Thus, we find here a combination of chemo- and thermo-diffusion contributions, which 
appears to be unique up to terms involving first-order space-time gradients. 



We observe that the diffusion current of Eq. (20) indeed generalizes the nonrelativistic 
ansatz discussed before, cf. Eqs. (0). For example, neglecting interactions, a straightforward 
calculation in the ultrarelativistic Boltzmann approximation gives: 

= rin[Cp,T-3] , (21) 

where C collects the constants. Using this in Eq. (pO]), i.e. fij = iJ,j{pj,T), we obtain: 

j; = -E-^.(^ + 3^) • (22) 

Obviously the term oc Vpj corresponds to the nonrelativistic one, while our covariant con- 
siderations additionally yield an unique thermo-diffusion term. 



Td^s^ = -Y,|la^-TY,(yr^{^^JZ^)^l{^,^Zj)^' 



Collecting the results of this section, we can now rewrite Eq. (14) with all the sources 
of entropy production in an explicitly nonnegative form: 

-^Q,A^Q,A^^ + |a^>V^A^A^W;:A^ + ad.un^ , (23) 

which is obtained with the help of Eqs. (p^-(p!S|) and by making use of the properties of the 
projector A. All transport coefficients are treated as constants here, despite the fact that 
they implicitly depend on the thermodynamical variables, if calculated microscopically pQ]. 



Together with condition (15), the Eq. (|23| ) guarantees that a closed system evolves from 
local towards global equilibrium. The latter is reached, of course, when all space-time 
gradients of T, fii, and u^^ vanish and the particle densities have their asymptotic values, 
such that the sources J7i vanish. 

Driving Towards Chemical Equilibrium. We now turn to the discussion of the local 
source terms J'i[pj], which describe the particle production in the continuity equations @. 
This will clarify the meaning of the constraint (p^). 

We assume a sufficiently small fluid cell with a given volume, such that all spatial gra- 
dients can be neglected in Eq. (^). This does not limit our considerations which are related 
to the invariant source terms J7i. Thus, in the local rest frame, we obtain approximately: 

dtPi = Ji = Y^[Gi^jPj - Lj^ipi] , (24) 
j 

where the source term is expressed in terms of gain and loss rates, G and L, respectively. In 
general, they are complicated function(al)s of the temperature and the densities themselves 



(see, for example, Refs. P, 15, M and Chap. X. in Ref. M). 



The master equations ( P4|) form a first order nonlinear flow system |23]. In the following 
we will show that Eq. (^) presents a sufficient condition for an attractive fixed point of the 
flow, which corresponds to the local chemical equilibrium. 

In order to proceed, we introduce the equilibrium densities Pi{T) = pi{fli,T), where 
p,i denotes the corresponding equilibrium chemical potential, and the reduced densities 
= Pi/ Pi, < (,i < oo. We assume that 

Pi{Ci,T)=il, + MCuT) , (25) 

with /j(l,T) = 0, is a strictly increasing function of at a given temperature T. This is 
exemplified by the result obtained in the Boltzmann approximation, Eq. (pl|). 



Pi{Ci,T)=fli + Tln[^i] . (26) 

In any case, the point = {(,i = 1, for all i} in the multidimensional ^-space corresponds 
to the chemical equilibrium state. 



Sufficiently close to equilibrium, ^ ^ ^i, the constraint (15) can be expanded. We 



suppress the temperature dependence from now on, tacitly assuming (as before) that local 
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thermal equilibrium is maintained by sufficiently fast elastic scattering among the con- 
stituents of the considered fluid cell. For convenience, we introduce rates weighted by the 
equilibrium densities, Gi^j = Git~.jPj and Lj^i = Lj^iPi. Then, employing Eq. (^) and 
using principle of detailed balance, 

E[G.-.(a)- V.(a)] = o , (27) 

i.e. i7i = in chemical equilibrium, we obtain: 

0<-^PiJi--Y. m- ii) ■ V{G,^j^, - Ij^.ii) . (28) 

i id 

The first nonvanishing contribution here is a positive gradient term; in the absence of 
equilibrium chemical potentials one obtains a curvature term here instead. Therefore, the 
nonnegative contribution —J2iPi>Ji to the entropy production, Eq. (p^), vanishes in the 
chemical equilibrium state. 

In order to prove that the equilibrium point ^ = is actually an attractive fixed point, 
we suitably rewrite the flow system, Eqs. (p^), 

dtiCr - 1) = Et^-.^J - Lj^^C^] , (29) 

j 

and introduce a distance measure M, 

i 

since f- = dfi/d^i > 0, pi being a strictly increasing function. Then, we obtain: 

dtMilii) = J2{p^{Ci) - -pi)-p^^t{^i - 1) = < o , (31) 

i i 

which holds sufficiently close to the equilibrium state, where Pi{Ci) ~ Pi ~ /i(l)('?i ~ !)■ We 
also made use of Eqs. ([T5|), (p^), and (|29|). Thus, the system evolves towards the chemical 
equilibrium fixed point, provided the constraint ( [l5| ) is satisfied. 

In order to arrive at the second equality in Eq. (^T]), we had to take into account charge 
and/or particle number conservation. The equilibrium values of the chemical potentials 
Pi are given by the equilibrium value of the chemical potentials associated to conserved 
quantities. In the case of a single charge conservation law this is written as: 

Pi = QiP , (32) 
9 



where qi denotes the generic charge of (anti-)particles of species "i". Thus we obtain the 
relation: 

fiiPidt(,i = fldt QiPi = , (33) 



which we employed in Eq. (31). The generalization including several distinct abelian conser- 



vation laws is obvious, while the consideration of nonabelian charges in this context presents 



an open problem, cf. Ref. [24|. 



In passing we mention that 

T-'^-dMMii,Ci)] , (34) 



may provide a useful measure of the chemical equilibration rate. According to Eq. (|31|) , it 
can be written as a sum of partial equilibration rates, = J^i'^^^^ with r-""*^ oc —piJi, 
while a more explicit form of the source terms, cf. Eq. (p^), allows to relate it to the 
equilibration times of chemical rate equations. 

Discussion of the First-order Formalism. Some general remarks about the first-order 
dissipative relativistic fluid dynamics presented in this work are in order here. It is labeled 
"first order", since the entropy current, Eq. (O) together with Eqs. (16)-(19), contains only 



up to first order gradients, which describe deviations from the equilibrium state of the fluid. 

It has been argued before that in a wide class of first-order theories the equilibrium 
state is unstable under physically admissible perturbations [25|. This includes the theories 



of Eckart [pq ] and of Landau and Lifshitz |27] especially. Unrealistically large growth rates 



of the unstable modes^ led the authors of Refs. to discard these theories in favor of 

the Israel- Stewart second-order theory [^], where only stable or damped modes are found. 

The expense is the increased complexity of the theory based on a larger number of degrees 

of freedom. This goes together with an increased number of phenomenological parameters, 

which ultimately should be calculated microscopically, as we discussed. 

We do not know, whether the first-order approach presented here prevents such unphys- 

ical instabilities. A full linear response analysis studying this problem is beyond the scope 

of this paper and is left for future work. 

^However, these were calculated only for macroscopic fluids, such as water. It is not obvious that the 
instabilities would be unrealistically fast for fluid parameters representing microscopic fluids, such as QGP 
droplets. 
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However, the growth of potentially unstable modes will be cut off by the finite size of 
the system in most realistic applications and by the additional dissipative processes taking 
place at the surface, such as radiation losses. These can be naturally incorporated in the 
future. Furthermore, we do take into account that the fluid may undergo dynamical changes 
of its composition. This is particularly important for relativistic fluids at extreme energy 
densities, where particle production is an important dissipative process, as long as chemical 
equilibrium is not fully established. Our considerations related to particle sources have 
shown that small perturbations of the chemical equilibrium are damped and decay. The 
system is driven by chemical reactions or particle production towards the equilibrium state. 

Conclusions. We have presented here the theoretical considerations which allow to deter- 
mine the entropy production which is consistent with the restriction to first-order space-time 
gradients of the thermodynamical and fluid variables, as well as with the requirements of 
the First and Second Laws of Thermodynamics. Our result in Eq. (|^) and the subsequent 
hydrodynamical framework are generally applicable to special relativistic fluids. As new 
elements, which need to be considered in applications to high-energy nuclear collisions and 
QGP, we introduce here the entropy producing mechanisms of chemo- and thermo-diffusion 
together with the contribution of particle production. 

We point out that nonlinear flow systems, such as the one describing the particle pro- 
duction processes in QGP, which we studied in relation to chemical equilibration, generally 
tend to show more structure than a simple attractive fixed point Therefore, one may 
speculate, whether some particular fluctuations may be caused by these nonlinearities in 
high-energy collisions (high multiplicity events), i.e. by the nonequilibrium chemistry of 
the reactions, with interesting consequences for the entropy production. 

Summarizing, the Eqs. (|l|), (|3|), and (^), together with the relevant definitions and sup- 
plemented by an equation of state, present a consistent framework for the study of dissipa- 
tive relativistic hydrodynamical systems. This comprises imperfect fluids with composition 
changing processes, particle production in particular, and the more familiar dissipative ef- 
fects of heat conduction, shear and bulk viscosities, and diffusion. We discussed the potential 
shortcomings of this formalism, which involves only up to first-order space-time gradients. 
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